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^ ' Abstract 

O-i. 

I \ We reanalyze the rare decays — > vr+z/i^ and J^l tt^vu in a general MSSM 
Qh. 

D . with conserved R-parity. Working in the mass eigenstate basis and performing adap- 



X 



tive scanning of a large space of supersymmetric parameters, 16 parameters in the 
constrained scan and 63 in the extended scan, we find that large departures from the 



5h ' Standard Model expectations are possible while satisfying all existing constraints. 

Both branching ratios can be as large as a few times 10"^'^ with Br{Ki^ -n^vv) 
often larger than Br{K~^ — > ir^i/v') and close to its model independent upper bound. 
We give examples of supersymmetric parameters for which large departures from the 
SM expectations can be found and emphasize that the present 90% C.L. experimen- 
tal upper bound on Br{K^ — > tt'^vu) gives a non trivial constraint on the MSSM 
parameter space. Unlike previous analyses, we find that chargino box diagrams can 
give, already for moderately light charged sleptons, a significant contribution. As a 
byproduct we find that the ranges for the angles jS and 7 in the unitarity triangle 
are relaxed due to the presence of new CP-violating phases in — and — B^ 
mixing to 12° < /3 < 27° and 20° < 7 < 110°. 



1 Introduction 



The rare decays — * ti^vv and Ki^ tt^i/u belong to the theoreticaUy cleanest pro- 
cesses in the field of meson decays. In fact their branching ratios can be computed to 
an exceptionally high degree of precision that is not matched by any other decay of 
mesons [1, 2, 3, 4, 5]. A detailed review of K ^ tiuu decays within the Standard Model 
(SM) and its extensions has been recently presented in [6], where an extensive hst of 
relevant papers can be found. Earlier reviews have been given in [7, 8]. 

As emphasized in [6], the clean theoretical character of these decays remains valid in 
essentially all extensions of the SM, whereas this is generally not the case for non-leptonic 
two-body B decays used to determine the CKM parameters through CP asymmetries 
and/or other strategies. While several mixing induced CP asymmetries in non-leptonic 
B decays within the SM are essentially free from hadronic uncertainties [9], as the latter 
cancel out due to the dominance of a single CKM amplitude, this is often not the case in 
extensions of the SM in which the amplitudes receive new contributions with different weak 
phases implying no cancellation of hadronic uncertainties in the relevant observables [10, 
11]. 

In this context an important virtue of — ttuu decays is the possibility of parametriz- 
ing the new physics contributions to their branching ratios, in essentially all extensions 
of the SM, in a model- independent manner by just two parameters [12], the magnitude 
of the short distance function X and its complex phase: 

X = |X|e^^^ (1.1) 

with \X\ = X{xt) and 6'x = in the SM. Here Xt = m^/M^. 

The most recent predictions for the relevant branching ratios within the SM read [6] 

Br{K+ 7r+i/z^)sM = (7.8 ± 1.2) • 10-^\ (1.2) 

Br(Ki^ 7r°i/P)sM = (3.0 ± 0.6) • 10'^^ (1.3) 

in the ballpark of other estimates [8, 13, 14, 15, 16]. As discussed in [6] a NNLO calculation 
of the charm contribution to K'^ — > tt'^i/P and further progress on the determination of the 
CKM parameters coming in the next few years dominantly from BaBar, Belle, Tevatron 
and later from LHC and BTeV, should eventually allow predictions for Br{K~^ 'k^vv) 
and Br{Ki^ — > tt^uu) with uncertainties of at most ±5%. 

On the experimental side the two events observed by the AGS E787 collaboration at 
Brookhaven [17, 18] and the additional event observed by AGS E949 [19] imply 

Br{K^ ^ TT+t/i/) = {U.7tlf) ■ 10"''- (1.4) 

While the central value in (1.4) is about a factor of two higher than the SM value, the large 
experimental error precludes any claims for signals of new physics in the n^uu data. 

Further progress is expected in principle from AGS E949, from the efforts at Fermilab 
around the CKM experiment [20], the corresponding efforts at CERN around the NA48 
collaboration [21] and at JPARC in Japan [22]. Hopefully this will bring additional 50-100 
events in the next five years. 
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On the other hand the present experimental upper bound on Br{Ki^ — > t^^uu) from 
KTeV [23] reads 

Br{Ki, ttVp) < 5.9 x lO"^. (1.5) 

This is about four orders of magnitude above the SM expectation but a Kj^ n^uP 
experiment at KEK, E391a [24], which recently took data, should in its first stage improve 
the bound in (1.5) by three orders of magnitude. While this is insufficient to reach the SM 
level, a few events could be observed if Br{Ki^ — > tt'^uu) turned out to be by one order 
of magnitude larger due to new physics contributions. Further progress on this decay 
is expected from KOPIO [25] at Brookhaven, and from the second stage of the E391 
experiment at JPARC [22], that could in principle provide a few hundreds of — > tt'^uP 
SM events, which would be truly fantastic! 

In this context let us recall that a model-independent upper bound on Br{Ki^ Tr^vi/) 
following from isopin symmetry reads [26] 

Br{Ki^ n^uu) < 4.4 • Br{K+ n+uu). (1.6) 

With the data in (1.4), which imply [19] 

Br{K+ TT^iyjy) < 3.8 ■ 10"^° (90% C.L.), (1.7) 

one finds then (1.6) 

Br(Ki^ 7r°i/p) < 1.7 • 10"'' (90% C.L.), (1.8) 

still two orders of magnitude below the upper bound from the KTeV experiment in (1.5). 

Recently it has been pointed out in [15, 16] that the anomahes seen in the B — > nK 
data may imply a significantly enhanced value of |X| and a large complex phase 9x' 

|X| = 2.17 ±0.12, = -(86±12)°, (1.9) 

to be compared with X = 1.53 ± 0.04 and = in the SM. 
The implications of this possibility are spectacular: 

Br{K+ TT+i/P) = (7.5 ± 2.1) • 10'^^ (1.10) 

Br(Ki^ 7r°i/p) = (3.1 ± 1.0) • 10"^°. (1.11) 

Consequently, while Br{Ki^ n^uv) ^ {l/?>)Br{K'^ 7r+i/p) in the SM, it is sub- 
stantially larger than Br{K^ — > Tr+i/p) in this scenario. The huge enhancement of 
Br{Ki, — > Ti^uu) found in [15, 16] is mainly due to the large weak phase 9x and its 
negative sign. One also finds 

Br(K+ TT+uP) ' ^ ^ 

which means that Br{Ki^ — > tt^vP) is rather close to its model-independent upper bound 
in (1.6). 

A spectacular implication of these findings is a strong violation of the golden minimal 
flavour violation relation [27] 

(sin2/3)^^p = {sm2(3)^Ks, (1-13) 
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where — /3 is the phase of the Vtd couphng that could in general be violated by new physics 
contributions to X — > tti/P and the asymmetry in Bd — ^ ipKs- 
Indeed, one finds [15, 16] 

(sin2/?),,p = sin2(/3 - 9x) = -{O.QQtol), (1-14) 

in striking disagreement with (sin 2P)^Ks — 0.736 ± 0.049. Other interesting imphcations 
are discussed in [15, 16, 28, 29]. 

The analysis in [15, 16] was model independent within a scenario in which new physics 
manifested itself only through enhanced EW penguins with a new complex phase. This 
scenario was first considered in [12, 30, 31, 32] in the context of rare K decays and the 
ratio e'/e measuring direct CP violation in the neutral kaon system, and was generalized 
to rare B decays in [33]. In [15, 16] it has been further generalized to S — > ttK decays. 

The question that we would like to address in the present paper is whether the spec- 
tacular effects in K ^ ttuu system found in [15, 16] are still possible within the general 
MSSM with new sources of fiavour violation when all known experimental and theoretical 
constraints on the relevant parameters are taken into account. 

New fiavour violating interactions are present in the general MSSM because the 
sfermion mass matrices can be non-diagonal in the basis in which all quark-squark-neutral- 
gaugino vertices and quark and lepton mass matrices are flavour diagonal. Instead of 
diagonalizing sfermion mass matrices it is often convenient to consider their off-diagonal 
terms as new fiavour violating interactions. This so-called mass insertion approximation 
(MIA) [34] has been reviewed in the classic papers [35, 36, 37], where further references 
can be found. 

In the context of the K — > ttz/z/ decays the most extensive analyses using the mass 
insertion method can be found in [38, 12, 30, 32]. It turns out, that sizeable enhancements 
of if — > TTi/P rates can only be generated by chargino-mediated diagrams with a large LR 
mixing in the up-squark sector. Reference [32], the most detailed of these papers, finds 
the upper bounds 

Br{K-^ TT+i/p) < 1.7 • 10-^°, Br{Kj^ ttVp) < 1.2 • 10"^°. (1.15) 

Larger values were found rather unlikely. Moreover, as discussed in detail in [12], in 
these models the minimal flavour violation relation in (1.13) can be violated due to the 
presence of a new phase 9x- A rough estimate shows that this phase could be as large as 
±25°. This is not as large as found in [15, 16] but sizable departures from the SM as seen 
in (1.15) can be found. 

However, in [12, 32] the following assumptions (approximations) have been made: 

• MIA (up to the second order of expansion) was used in calculation of if — > nuP 
decay rates and other related processes. 

• limited number of MSSM parameters, assumed to be most important, were used in 
numerical scans. 

• bounds on mass insertions coming from various experimental results were obtained 
by the approximate procedure of requiring that each individual term in the mass 
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insertion expansion at most saturates the measured value; this neglects the possi- 
bility of significant cancellations between contributions from different terms in the 
expansion and also interference with the SM contribution that is always present. 

• some additional assumptions not coming from low-energy physics, like GUT induced 
RGE relations between mass insertions, were taken into account. 

The questions that we would like to address here are then, whether by relaxing the 
assumptions made in [12, 32], going beyond the MIA and performing more extensive 
numerical scans, while satisfying all existing experimental constraints on supersymmetric 
parameters, 

• 9x can be as large as found in [15, 16], 

• Br{K^ Ti^vv) can be significantly enhanced over the SM expectations so that it 
is at least as high as its central experimental value in (1.4), 

• Br{Ki^ — > 7r°i/p) can be enhanced by an order of magnitude with the ratio Br{Ki^ — > 
7r'^uP)/Br{K'^ — > n'^uP) reaching the bound in (1.6). 

Moreover, there are the questions: 

• In what regions of the parameter space of the general MSSM do the maximal de- 
partures from the SM expectations occur? 

• What are the bounds on the MSSM parameters implied already by the experimental 
value on Br{K^ — >■ ti^ui') in (1.4)? 

Answering these questions is a non-trivial numerical task, due to the large number 
of free parameters and experimental constraints which have to be considered. Here we 
would like to demonstrate an efficient method of a random scan over the MSSM parameter 
space, based on an adaptation of the Monte Carlo integration algorithm VEGAS [39, 40]. 
Such a method is designed to automatically concentrate most of the randomly generated 
points in the SUSY parameter ranges giving the largest deviations from the SM results, 
thus allowing for analyzing very large parameter spaces, with 20 or more dimensions, in 
a reasonable time and without very extensive computer CPU usage. This allows us to 
answer the first three of the questions listed above. The remaining two points will be 
discussed elsewhere. 

Our paper is organized as follows. In Section 2 we will present the formulae for the 
branching ratios Br{K^ n^uu) and Br{Ki^ -k^vv) in terms of the function X, 
discuss their properties and compare the general MSSM considered here with the simple 
new physics scenario analyzed in [15, 16]. In Section 3 we discuss the problematic of the 
huge space of parameters in supersymmetric theories and the method that we propose to 
tackle this problem. In Section 4 we present our numerical results. Finally, in Section 5 
we summarize the main findings of our paper. A list of the contributions to the function 
X from the relevant box and penguin diagrams is given in an Appendix. 
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2 Basic Formulae 



2.1 Effective Hamiltonian 

The effective Hamiltonian relevant for n^uu and Ki^ tt^vi' decays in a general 

MSSM considered in our paper can be written as follows 

where the internal charm part 

^2= E V:,V,A{sd)v-A{m)v-A (2.2) 

l=e,iJ,,T 

is fully dominated by the SM contributions and 

^2= E VtlVtd[XL(sd)v-A('^ii^i)v-A + Xn{M)v+A(i^ii^i)v-A] (2.3) 

l=e,iJ,,T 

with Xl receiving both the SM and supersymmetric contributions and only the latter 
ones that also include charged Higgs exchanges. In the SM X^ is real and given for 
rnt{mt) = (168.1 ± 4.1) GeV by 

Xl^ = XsM = 1.53 ± 0.04 . (2.4) 

The index I = e, n, t denotes the lepton flavour. The dependence on the charged 
lepton mass resulting from the box diagrams is negligible for the top contribution. In the 
charm sector this is the case only for the electron and the muon but not for the r-lepton. 
Below we will give the branching ratios that follow from (2.3). 

This parametrization of is very useful for phenomenological applications and for 
the comparison with [15, 16], but one should remember that in a general MSSM scenario 
considered here not all contributions to arc proportional to V^*Vtd. This means that 
parametrizing Tf^^ as in (2.3), necessarily puts some CKM dependence into Xl and Xr. 
Therefore we will use below also a different parametrization of new physics contributions. 

Now, as the strong interactions are not sensitive to the sign of 75, the hadronic matrix 
elements of {sd)v-A and {sd)v+A are equal to each other, and consequently the function 
X playing the central role in our paper is simply given by 

X ^ Xl + Xr^ |X|e^^^ . (2.5) 

Explicit expressions for Xl and Xr are collected in the Appendix. As one can see there, 
chargino and part of the neutralino contributions to Xr are strongly suppressed by small 
Yukawa couplings of the down quarks. As we will see, the remaining neutralino con- 
tribution proportional to U{1) gauge coupling is typically also much smaller than the 
dominant terms in X^, so that \Xr\ <^ [X^l also in the MSSM (except may be for some 
non-interesting points where X^ is small due to cancellations) and can be neglected for 
all practical purposes. 
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2.2 Branching Ratios 

The branching ratios for the K 
can be written as follows [6] 



Truiy decays, in the new physics scenario considered, 



A5 



A 



A5 



/im(AiX) 



A5 



(2.6) 



(2.7) 



with A = 0.224 being one of the Wolfenstein parameters [41], A^ — V^gVtd, Ac — V*gVcd, 
K-factors equal to [6] 



= (4.84 ± 0.06) • 10 



-11 



kl = (2.12 ± 0.03) • 10 



-10 



and 



2^NL + 3^NL 



0.39 ±0.07 



(2.8) 



(2.9) 



resulting from the NLO calculations in [4, 5]. The anatomy of the error in Pc{X) has 
been recently presented in [6]. 

As discussed in the next section we will use as our input parameters 



IK.I = A = 0.224, 



\Va\ = AX'^ = 0.0415, Rb = 0.37 



where 



Rh 



:i - AV2) 

A 



ub 



cb 



(2.10) 



(2.11) 



We recall that A is a. Wolfenstein parameter [41] with A = 0.83 ± 0.02 and is one 
of the sides of the unitarity triangle. 

While the parameters in (2.10) contain uncertainties, the latter are sufficiently small 
so that they can be neglected in comparison with hudgc uncertainties in the values of 
super symmetric parameters. The choice of these three parameters is dictated by the 
fact that they are extracted from tree-level decays and consequently their values are not 
subject to new physics uncertainties. 

As the fourth variable we will take the angle 7 in the unitarity triangle that is equal to 
the CKM phase Sckm in the standard parametrization of this matrix. The angle 7 can be 
in principle measured without any new physics pollution in tree-level B decay strategies 
that will be only available at LHC and BTeV [9]. In the general MSSM the value of 7 
may deviate from the one extracted from the usual analysis of the unitarity triangle that 
uses SM expressions. For this reason we will allow 7 to vary in the full range 



■180° < 7 < 180°, 



(2.12) 



but as we will see in Section 4, only the range 20° < 7 < 110° is allowed when all 

constraints are taken into account. 
We define next Xsusy through 



(2.13) 
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and introduce 



|2 



a = PciX) + ^-^XsM ~ 1.43, b = rJ-^Xsm ~ 0.39, (2.14) 
that do not depend on supersymmetric parameters. We find then 

Br{K+ ^ TT+i/P) = K+ [{Qeffil) + ImXsusYf + {PefAl) - ReXsusyf] , (2.15) 

Sr(XL ^ Ti^'yV) = «i (QefAl) + ImX5^/5y)' , (2.16) 

where 

PefAl) = (1 - y) [a - 6COS7] , Qe//(7) = &sin7 . (2.17) 

As Peffil) and Qeffil) can be fully determined within the SM, provided 7 can be mea- 
sured through tree-level decays, the formulae (2.15) and (2.16) transparently exhibit su- 
persymmetric contributions. With 7 65° we have 

Pe//(7) ~ 1-23, Qeffil) ^ 0.35 . (2.18) 
2.3 Comparison with [16] 

In order to compare our results with the ones in [16], it is useful to have the expressions 
for both branching ratios in terms of Rt and /3, where (3 is another angle and Rt another 
side of the unitarity triangle. 

Now, for each pair (i?6,7), used above, one can determine the pair {Rt,P) through 

Rt = Jl + Rl-2RhCos^, cot/3 - ^'^bcos^ ,2.19) 
* Rb sin 7 

and consequently calculate both branching ratios by using 

Br{K+ TT+i/i?) = K+ [f^A'^RflXf + 2fPciX)A'^Rt\X\ cosPx + Pc{Xf] (2.20) 
Br(Ki^ 7r%P) = KLr^A^R^\X\'^ siii^ (3x, (2.21) 

where 







(2.22) 




f = 


Vts 


« 0.985, 


(2.23) 


Vu = ARtX'e-'^, 


Vts = 






(2.24) 



The small phase j3s ~ —1° can be neglected for all practical purposes. 

While the formulae (2.20) and (2.21) look exactly as the ones given in [16], there 
is one important difference that should be empasized here. In the new physics scenario 
considered in [16] the values of /3, Rt and 7 were unaffected by new physics contributions 
and the sole effects of new physics were felt only through \X\ and 9x- In a general MSSM 
the a^Ks asymmetry does not generally measure (3 but rather (3 + 9d with coming from 
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new complex contributions in i?^ — i?^ mixing. Thus the constraint on the angle (3 coming 
from the a^^s asymmetry reads 

sin 2(/? + Od) = 0.736 ± 0.049 (2.25) 

so that 13 could be quite different from (3 — 23.5° found in the SM and used in [16]. 
Moreover Rt can be modified. Therefore, the comparison of the results obtained here 
with the ones in [16] requires some care. 



3 Numerical analysis 

3.1 Independent MSSM parameters 

In the general MSSM, the predictions for the branching ratios considered here and various 
experimental constraints can depend on almost every MSSM parameter. The dependence 
on the majority of them is weak enough, however, to be safely neglected in principle for 
the purpose of computing Br{KL — > tt^vP) and Br{K^ — > n^uV). 

Still, parameters may indirectly enter through the expressions for other observables, 
which we use to constrain the parameter space. Taking a conservative approach we 
allow for the independent variation of more free parameters than included in the existing 
literature. 

Kobayashi-Maskawa phase -180° < 7 < 180° 

CP-odd Higgs boson mass 150 < < 400 

SU{2) gaugino mass; we use Mi GUT-related to M2 50 < M2 < 800 

Gluino mass 195 <m~g< 2000 

Supersymmetric Higgs mixing parameter —400 < /i < 400 

Common fiavour diagonal slepton mass parameter 95 < Mgi < 1000 

Common mass parameter for the first two generations 

of squarks, as well as Ir 240 < Mgq < 1000 

Squark mass parameter for the left stop and sbottom 50 < M^^ < 1000 

Squark mass parameter for the right stop 50 < M^^ < 1000 

Flavour universal trilinear scalar mixing parameter 

(normalized to the fermion mass) — 1 < ^ < 1 

Mass insertion Sl\ \S^^l\ < . 1 35 

Mass insertion I^ijlrI < 1-65 

Mass insertion 5^^^ \^i?lr\ < 1-65 

Table 1: Parameters and their ranges used in the "constrained" scan (16 real degrees of 
freedom). All mass parameters are in GeV. 
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To start with, we assume that Br{KL — 7r°i/i/) and Br{K^ — > tt'^uu) may depend 
significantly on the set of unknown SM and SUSY parameters hsted in Table 1. We took 
as fixed the other SM parameters, including fermion masses and the CKM parameters 
given in (2.10), as their measurements are known to be relatively insensitive to physics 
beyond the SM. This leaves as our only free SM parameter the CKM angle 7, which to 
date has not been determined from tree-lcvcl decays. 

Apart from the flavour diagonal super symmetric mass parameters, we vary several 
off-diagonal ones. For a transparent parameterization and easy comparison with the lit- 
erature, it is useful to define mass insertions, even though our actual computation makes 
use of exactly diagonalized sfermion, chargino, and neutralino mass matrices. Follow- 
ing [35, 36], we define 

Here /, J denote quark flavours, X, Y the (supcrflcld) chirality, and the "sector" S can 
be either U or D. Evidently, this definition depends on the basis chosen in fiavour space. 
For the right-handed up and down squarks, there exists a preferred choice selected by the 
mass eigenstates of their respective SM partners. In the left-handed sector this is different, 
because SU (2) invariance dictates that the MSSM Lagrangian contains only one soft mass 
matrix for the doublet squarks. This implies that ^Ifj^]^ and are related. In the super- 
CKM basis, where each squark field is the superpartner of a quark mass eigenstate, this 
correspondence is given by CKM matrix elements. As a consequence, even for diagonal 
{^dYll off-diagonal 5^^^ are in general nonzero in this basis. Also, the mixings and 
mass splittings among left-handed up squarks are constrained by bounds on down-squark 
mixings from — mixing. 

In the context oi K ^ ttuu, it has been pointed out by the authors of [12] that a 
particularly useful choice of basis is one where the left-handed up squark basis states 
are simply the SU{2) partners of the left-handed down squarks, which in turn are the 
supcrpartners of the left-handed down-type quarks, so that the gluino and neutralino 
couplings are fiavour conserving. This eliminates the CKM elements from the chargino 
couplings to the left-handed up squarks, and in particular from the term in the MIA 
that is believed to be dominant [30] . Besides, such a basis ensures that any choice of 5's 
automatically preserves SU (2) invariance. We use this basis, to be called the BRS basis, 
in presenting our results in this paper. We also omit the "sector" index on the Sll- 

The decays Kl ti^vv and n^uu arc AS = 1 processes and correspond to 

a transition between the first and second quark generations. Thus, they could certainly 
depend on (12) squark mass insertions. However, as pointed out in [30], both decays are 
sensitive also to second-order terms in the mass insertion expansion, namely, to products 
of a (13) and (the conjugate of) a (23) mass insertion. Therefore, we start from varying 
independently five (12) mass insertions, ^"dlr ^"^^ ^ij'lr^ '^\v,s five more 

(13) and five (23) mass insertions (assuming for the moment that all LR mass insertions 
are hermitian, e.g. 6]^lr = ^dlr)- 

Even assuming that all flavour diagonal quantities from the list above are real (there 
is no reason to constrain mass insertions to be real), we have listed already 40 free real 
parameters. This is a huge number for any reasonably dense numerical scan. To avoid 
excessive computation time, we flrst tested how sensitive to the various mass insertions 
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the predictions for our branching ratios really are. In Fig. 1 we plot the dependence of 
Br{K'^ — > TT+i/f?) on 5}^]^ and the product 5}j]^ji5'ijlji for a chosen set of SUSY parameters. 
The dependence on the other mass insertions, not shown in the plots, is much weaker and 
can be neglected in the first approximation. 



X 10 



X 10 
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IOg(6uLR(5uLR)/2 



Figure 1: Dependence of Br{K^ ti^uu) on chosen mass insertions for 7 = 65°, tan (3 = 
2, Msg = Mi^ = Mi^ = 500, Msi = 300, M2 = 200, = 3M2, fi = 100, Ma = 200 
(all masses in GeV). Solid, dashed and dotted lines: and ^[/LR^^Li^ complex, real and 
imaginary, respectively. 



Fig. 1 suggests that it is sufficient in the numerical scan to vary only three independent 
mass insertions, 5f7^2._R 5^^/?. In fact the dependence on which is tightly 
constrained from mixing in the neutral kaon system, is almost negligible for realistic 
values. 

Also, as we tested, both branching ratios in question are only weakly sensitive to 
GUT- or SUGRA-type assumptions concerning relations between SUSY mass parameters 
(e.g unified gaugino masses or sfermion mass parameters). Thus, the large but already 
numerically more feasible total number of 16 free parameters appears to be sufficient to 
explore the possible ranges of Br{K^ — > tt^vD) and Br{KL tt^vV) in the framework 
of the MSSM. However, this is not true when trying to satisfy experimental bounds and 
selecting allowed values of the SUSY parameters. For example, keeping the product 
^^ll^Ill constant but varying independently does not affect Br{K~^ tt+z/z/) and 
Br{KL TT^i/i/) but allows to satisfy bounds coming from AMd, the measured B^ — B^ 
mass difference. This motivates keeping, in Sj^^, one more flavour violating parameter 
even in the "minimal" 16-parameter scan. Later, we will see that our numerical method 
is powerful enough to study the consequences of further increasing the number of degrees 
of freedom, which turn out to be unimportant. 
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3.2 Theoretical and experimental bounds 

The enormous freedom in the SUSY parameter space is reduced by a number of phe- 
nomenological constraints the theory must satisfy. From the experimental side we take 
into account the set of bounds and measurements hsted in Table 2. As can be seen, we do 
not take into account the e' je constraint, as its calculation still has very large theoretical 
uncertainities. Doing this we are aware of the fact that e' je could one day offer a very 
powerful constraint on the size of SUSY contributions [31, 32]. 



Quantity 


Measured value 


Experimental error 


Lightest neutralino mass 


> 46.0 GeV 




Second lightest neutralino mass 


> 62.4 GeV 




Lightest chargino mass 


> 94.0 GeV 




The two "sbottom" masses 


> 89.0 GeV 




The two "stop" masses 


> 95.7 GeV 




all other squark masses 


> 250.0 GeV 




\^k\ 


2.280 ■ 10-3 


0.013 ■ 10-3 


AMk 


3.489 • 10-15 GeV 


0.008 • 10-15 


AM, 


3.31 • 10-13 GeV 


0.04 • 10-13 


AM, 


> 9.5 • 10-12 GeV 




Br{B, ^ X,7) 


3.28 • 10-^ 


+0.41 -, r,-4 

-0.36 ■ 


(sin2/3)^Ks 


0.736 


0.049 



Table 2: Experimental measurements used to constrain the MSSM parameter space. Lim- 
its on supcrsymmctric particles masses published in [42] are used. The quotes on "sbot- 
tom" and "stop" indicate that these are only approximate flavour states. 

Due to the large number of unknown parameters and substantial theoretical uncer- 
tainities present in many calculations, it is pointless (and certainly very difficult to im- 
plement from a practical point of view) to use everywhere the state-of-art calculations of 
higher order QCD corrections and other non-leading effects calculated in the literature. 
Nevertheless, whenever feasible we try to avoid unnecessary simplifications made in many 
papers. In particular: 

• For low energy FCNC and CP violating processes, we do not restrict ourselves to 
the most commonly used MIA. Instead, we calculate all Wilson coefficients relevant 
for a given process in the mass cigcnstate approach, taking into account the full set 
of contributions - SM, Higgs, chargino, neutralino and gluino sectors. Further, we 
compare experimental results with the full 1-loop SUSY expressions, not just the 
dominant (usually gluino) term. 

• To complete the bounds coming from the low-energy data, we also compare lower 
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mass bounds for SUSY particles, obtained from accelerator experiments, with the 
tree-level eigenvalues of the MSSM mass matrices. 

For comparing theoretical predictions for the low-energy observables with experiment, 
we apply the following procedure. For a given set of SUSY parameters, we calculate 
all appropriate parton-level diagrams. Next, we construct expressions for the considered 
quantities using the central values of the QCD evolution factors and necessary hadronic 
matrix elements, obtained by perturbative SM and lattice-QCD computations, respec- 
tively. Finally, for every quantity Q we require: 

igexp _ gtfe| < ^/^Qexp ^ ^|gt/.|^ ^ ^ g.S (3.2) 

with the exception of AM^, for which we require (1 + > Q'^^^. 

The first term on the RHS of (3.2) represents the 3a experimental error. The second 
term corresponds to the theoretical error. In principle, such an error differs from quantity 
to quantity and is usually smaller than q = 50%, which we assumed as a generic number 
in all calculations. However, apart from the theoretical errors coming from uncertainities 
in the QCD evolution and hadronic matrix elements calculations, one should take into 
account also problems arising due to the limited numerical scan density. In principle, with 
a very dense scan, it should be possible to find SUSY parameters fulfilling (3.2) within the 
"true" theoretical errors of present calculations. Such a dense scan requires, however, a 
huge amount of computer time - with 16 or more free parameters and rather complicated 
mass eigenstate formulae, it would take months of CPU time. This does not seem to be 
necessary and may even be undesirable. Our goal is to find "generic" allowed values for 
the K — > nuP decay rates, i.e. values possible to obtain for fairly wide ranges of SUSY 
parameters, without strong fine-tuning and resorting to some very particular points of the 
parameter space where the experimental bounds are satisfied due to precise cancellations 
of various types of contributions. Thus, in our scan we use wide "theoretical" errors, 
assuming that this procedure points to the correct ranges of the MSSM parameters, and 
if necessary the exact values of such parameters fulfilling the bound in (3.2) with smaller 
q could always (or at least almost always) be found. This could be achieved e.g. by a 
denser scan or by more advanced numerical routines, possibly even solving numerically 
the set of non-linear equations given by all constraints. 

The only exception from the procedure described above is the imposition of the con- 
straint from the CP asymmetry a^i^g(t). Making the safe assumption that the supersym- 
metric contributions to the decay amplitude can be neglected, this asymmetry measures 
the phase of the — amphtude 

= (S°|i/e//(A5 = 2)15°) = |Mi2|e^^(^+^'') (3.3) 

through 

ai,Ks (t) = -a^Ks sm(AMdt) = sin 2{(3 + Oa) sin( AM^i) . (3.4) 

Here — /3 is the phase of Vtdi while Od is the new effective phase coming from the super- 
symmetric contributions. As in the general MSSM the presence of non-vanishing sfermion 
mass insertions implies contributions of several local operators to M12, the phase 9^ suffers 
from potential uncertainities related to the hadronic matrix elements of these operators. 
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Only if one of these operators dominates over the others in the full amplitude can 9d be 
cleanly related to supersymmetric parameters. In other cases its error can be sometimes 
very large and moreover is hard to estimate - it would require additional scanning over 
hadronic uncertainities for each given set of SUSY parameters. To avoid such problems, 
in our analysis we assume conservatively 



exp 



sin2(/3 + ^,)|<2Aa;^4 



(3.5) 



with sin2(/3 + 64) calculated in the MSSM for central values of the hadronic parameters 
0.736 ± 0.049. Expressing everything in degrees we require then 



and 



|23.7° - (/3 + ^d)| < 4.2°. 



(3.6) 



Assuming only experimental and no theoretical errors on the R.H.S. of (3.5), we reject 
some otherwise valid points from our scan, but we checked that this does not have any 
significant effect on the results discussed in following sections. 

Apart form the experimental bounds, there arc also bounds from the requirement that 
the vacuum is stable, or that the true ground state of the theory does not break color and 
charge. We apply the corresponding Charge and Color Breaking (CCB) and unbounded 
from below MSSM scalar potential (UFB) bounds [43] , which give constraints on the A- 
terms and consequently on the left-right elements of the sfermion mass matrices. We use 
the tree-level expressions in the form: 



AlJ 
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< 
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< 



< (1 + 25^^)/|y7-(^.-^)p + {miy^ + + |// 
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+ {ml] 


KK' 







(3.8) 

As suggested in some papers (see e.g. [44] for a comprehensive discussion), tree-level 
CCB/UFB bounds may not be very stable against radiative corrections. Again, we do 
not try to use more elaborate formulae, but relax the tree-level bounds by adding factors 
/ = 2 on the RHS of all equations in (3.7) and (3.8). (Strict tree-level bounds correspond 
to / = 1.) 



3.3 The adaptive scan algorithm 

The numerical analysis laid out in the previous subsections requires us to scan an A^- 
dimensional (A^ > 16) parameter space sufficiently densely to account for all regions 
where the quantity X of (2.5) possibly is large. A straightforward approach is to scan 
over a uniform A^-dimensional grid. However, this is both time-consuming and ineffective. 
Nevertheless, we have performed some test scans of this kind to get a first notion of the 
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attainable rates of both branching ratios considered, as well as of X. We generated about 
1.6M points, only about 1% of which passed all theoretical and experimental constraints. 
The distribution of the remaining points is shown in Fig. 2. What can be observed is 
that most points are close to the SM value of X, i.e. new-physics effects are very small 
for most grid points. The departure from the SM lies along straight lines, which is an 
obvious artifact of the small grid density, in spite of the fact that more than 1 million 
points were generated. 

However, what we are interested in is not so much the somewhat vague notion of 
a "probability density" for quantities shown in various scatter plots, which is clearly a 
function of the choice of parameters, of their scaling chosen during the scan (e.g. linear 
or logarithmic) and the range considered for each parameter. Rather, we are interested 
in the maximal allowed regions of the parameter space and the extremal values of the 
observables, viz. the (complex) function X, which are physical (i.e. parameterization- 
invariant). To find this information, a dense and very time-consuming scan is necessary. 



***** * 



0.02 0.04 0.06 0.08 



0.12 0.14 0.16 
Br(K*->7vVi/) 



Figure 2: Scatter plot of X, Br(K^ n^uu) and Br^Ki tt^uu) distributions for 
tan/9 = 2 and uniform grid scan. Parameters varied in the ranges < 7 < 180°, 200 < 
< 500, 150 < Msi < 300, -400 < fi < 400, 200 < M2 < 600, 150 < Ma < 300, 
— 1 < ^4 < 1, mass insertions between 0.003 and 0.3, their phases set to 0°, 45° or 90°. 



In practice, due to the number of free parameters, the grid involves only 2 to 5 different 
values for each parameter, posing an additional risk that a relevant region of the parameter 
space may fall through the grid altogether. In fact, this will be likely if the largest effects 
do not simply occur where some parameters are smallest and others are largest, but in 
areas where some type of cancellation, even moderate, is needed to fulfill the experimental 
constraints. 

A random generation of parameter points may exhibit some of these overlooked regions. 
Such an approach is particularly useful when, as it is true for the quantities considered in 
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this paper, the analyzed expressions depend strongly only on a subset of M free parame- 
ters, M < N. In this case one can think in terms of an "effective density" - the projection 
of the random points onto the M-dimensional subspace of "relevant" parameters now has 
a density (in that subspace) that is about (2 . . . 5)^/*^ times higher than for a uniform 
grid. To see this, envision projecting the original grid: for given relevant parameters, the 
remaining N — M parameters describing a given point, which are approximately irrele- 
vant, just parameterize points projected onto the same M-dimensional point. Thus in the 
uniform-grid approach one repeats the same analysis many times. Conversely, the projec- 
tion of the randomly distributed points is again randomly distributed over the subspace, 
and each new point gives additional information. 

Note that in general the subset of relevant parameters is different for each observable. 
This precludes just removing the irrelevant parameters from the scan, in the case one be- 
lieves to know them beforehand. The random scan, on the other hand, gives an improved 
effective density for each of the corresponding subspaces. 

Still, for randomly distributed parameter points, an expression \X — X'^^^ describing 
the deviation from the SM prediction (or similar expressions for other observables, which 
we use to constrain the parameter space) will give a small number over most of the 
parameter space. On the other hand, in the regions where it is large it also tends to vary 
more strongly. Therefore, while effectively improving the grid density, random generation 
of points with a uniform distribution may still miss the exact maxima for realistic numbers 
of points. It would be desirable to generate points more densely where the function X 
is large and/or varies quickly - such a procedure would be particularly efficient if, as is 
expected, our results mainly depend on just a few parameters. 

Recently, Brein [39] has suggested using an adaptive Monte-Carlo integration algo- 
rithm such as VEGAS [40] for a similar purpose. VEGAS performs "importance sam- 
pling" via an iterative algorithm. At each iteration, it generates a certain number of points 
according to a probability distribution determined by the integrand values encountered 
during the previous iterations. The probability distribution is chosen to be separable and 
is adjusted after each iteration. The initial distribution is chosen to be uniform. This 
procedure is designed to minimize the statistical error of the integration by increasing the 
number of points in those parts of the integration volume where the integrand is found to 
be large. While we are not interested in computing any integral to any given precision, 
VEGAS does provide what we desire, if we choose our integration volume as the parame- 
ter space (suitably rescaled) and the integrand such that it becomes large for large values 
of X. We found the following "integrand" useful: 

{0 parameter set rejected by constraints 

X-X^^P constraints satisfied ^^'^^ 

where X is the quantity defined in (2.5), X^^ is its SM central value (2.4) and we varied 
the power n between 3 and 8. (We do not write explicitly the dependence of X on 7 

and SUSY parameters.) We "integrated" this function numerically with the VEGAS 
routine, storing all generated Monte Carlo points in a separate file along with the values 
of various observables for analysis. The parameter ranges used for the integration are 
given in Table 1. 

Using the VEGAS algorithm has the great advantage of sampling mostly the important 
regions of parameter space, where the function X really depends significantly on at least 
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some of the parameters. It also allows us to increase safely the number of the degrees of 
freedom — adding new variables has only a moderate effect on total computation time as 
long as the function X, or in general the VEGAS integrand, is weakly dependent on them. 
In fact, we were able to perform Monte Carlo sampling over huge portions of the MSSM 
parameter space, probably never tried before, of up to 63 dimensions (parameters), and 
to judge only afterwards, from the obtained distributions, which of the parameters where 
important for a given problem. 

The applied procedure carries some dangers. First, one should mention a point that 
is not always appreciated. As is well known, the "random" numbers provided by various 
standard libraries are in fact deterministic sequences which may or may not satisfy cer- 
tain conditions of randomness, or absence of correlation. In our context, when one uses 
N subsequent numbers in the sequence to define a point in an A'^-dimensional parame- 
ter space, the widely used generators based on linear congruences tend to concentrate 
these points on a bundle of (A^ — l)-dimensional hyperplancs. This can reintroduce the 
danger of points of interest "falling through the grid" mentioned above. To address this 
problem, we used two different random number generators, both of which avoid such cor- 
relations. The first is a subtractive Mitchell-Moore generator discussed by Knuth [45] in 
an implementation of Kleiss [46] . The second is based on a combination of Fibonacci and 
arithmetic sequences proposed by [47] in an implementation of James [48]. We obtained 
similar plots for both algorithms, which strengthens our conviction that our procedure is 
stable and does not miss any relevant regions. 

Also, there is the question of stability of data samples obtained in this way. In order 
to check it, we varied the power n in (3.9) in the range of 3 to 8, and inspected how the 
shape of the plotted regions evolved with the numbers of points generated. The choice of 
the particular function / appears to affect mostly the time in which the shape of various 
distributions stabilizes, not the boundaries itself. Similar comments apply to tuning of 
the internal control parameters of the VEGAS routine. 

4 Results 

4.1 Sensitivity to scan parameters 

As a first step of our analysis, we plotted distributions ol\X\ obtained in the "constrained" 
scan (described in detail in Section 3.1) versus various flavour diagonal and off-diagonal 
parameters, in order to check the sensitivity of the considered decay rates to them. The 
plots obtained support our conjecture that only a modest subset of the MSSM parameters 
is in the first approximation relevant for our analysis, even after including the compli- 
cated set of additional bounds. Most of the obtained distributions for fiavour diagonal 
parameters are flat (dependence on M4, Afi, nig) or almost flat (dependence on Af^/, A). 
A more pronounced dependence could be observed only in the case of 7 and the SUSY 
parameters plotted in Fig. 3. 

The dependence visible for the other parameters can be explained by looking at Fig. 4, 
where we plotted the contributions from the chargino and neutralino sectors to \tX. 
The chargino diagrams are always strongly dominant, typically one order of magnitude 
larger than the neutralino contributions. Gluino exchanges (in our case penguin diagrams 
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X 




(GeV) Id. (GeV) 



Figure 3: Distributions of X versus 7, Mf^, M2 and ji for tan /3 = 2 in data set generated 
during VEGAS integration. 
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only) can be always completely neglected, as we checked that they are several orders of 
magnitude smaller than other contributions. Thus, the size of both decay rates, Br{K^ — > 
n'^yV) and Br{KL — > tt^i^P), is determined by the chargino-up-squark contribution and 
should depend mostly on the parameters entering the expression for this amplitude, in 
agreement with the results of [12, 30, 32]. This is not true for the variety of experimental 
and theoretical bounds we took into account in constraining the MSSM parameter space 
- they also depend on other parameters (e.g., the hole in the distribution comes of 
course from the bounds on the lightest neutralino and chargino masses). This has some 
secondary influence on the shape of the plots in Fig. 3, as the imposed constraints can 
lead to correlation between the allowed ranges of some of the parameters directly relevant 
for the K — > ■nvi' decay calculation. For instance, attainable values of |X| grow first with 
MSSM masses like M^^ or M2, because experimental constraints for squark mass insertions 
are easier to satisfy for heavier SUSY particles (similarly for distributions with Mi^ or 
Mgqi not included in Fig. 3). Later they go down again because the suppression of SUSY 
loop diagrams contributing to K ^ ttvi) decays for heavy virtual particles dominates 
above the effect of weakening the impact of the experimental constraints. One should 
also note that the distribution of 7 in Fig. 3 shows an increased density of points around 
SM-preferred value 7 ^ 70° ± 30° - for the remaining points, the larger or smaller value 
of the KM phase must be compensated by the phases of some mass insertions. 



X 

0.5 - 
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Figure 4: Chargino and neutralino contributions to A^X for tan/3 = 2. 



An interesting observation can be made by comparing the relative size of penguin- and 
box-type contributions to the K — > ttj/P decay amplitude. A common assumption used in 
the literature says that box diagrams are parameterically suppressed by 0{M^r/M^i^p^^^) 
and can be safely neglected when compared with Z'^-penguin contribution. This is cer- 
tainly not justified for light slcpton masses, and not very accurate even for moderate 
Msi ~ 300 GeV. In Fig. 5 we plot the absolute value of the ratio of box to penguin contri- 
bution against the slepton mass, including in the distribution only particularly interesting 
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points for which the — > t^^uu decay rate is large, Br{K^ — > t^^uu) > 1.5 ■ lO"^''. As 
can be immediately seen, for slepton masses just above 100 GeV box and penguin con- 
tributions are in a substantial number of cases comparable in amplitude, and even for 
Mgi = 300 GeV their ratio can still reach 30%. Thus box diagrams definitely should be 
taken into account in realistic calculations. Of course, their presence introduces slepton 
mass dependence into the considered branching ratios, which would otherwise be negligi- 
ble. 




Ms, (GeV) 

Figure 5: Absolute value of the ratio of box to penguin contributions as a function of 
slepton mass for tan/3 = 2. Only points for which Br{K^ — > t^^vV) > 1.5 • 10~^° are 
plotted. 

In our "constrained" scan we varied 3 types of mass insertions, 6}^j^, and We 
checked that the dependence on 5}^^ is actually almost negligible, as expected because it is 
strongly constrained by bounds coming from and AM^ measurements: < 0.05. In 
Fig. 6 we plot the dependence on the moduli of the remaining two ULR mass insertions. 
The dependence on them is quite pronounced, in agreement with the conclusions of [30] 
that the second order ULR terms in the MI expansion give the dominant contribution to 
the considered decays. We return to this point in Subsect. 4.3 below. For comparison, we 
also plot \X\ against just the real part of 5^^^,^, which exhibits less correlation. 

4.2 Allowed ranges for the X function, Br{Ki, — > tt^uu) and 

The most interesting information from scanning over SUSY parameters are the maximal 
values of \X\, Ox and the branching ratios Br{KL — > 7r°i/p), Br{K^ — > tt+i^P) which can 
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Figure 6: Distributions of X versus real part and modulus of 5ulr niass insertions for 
tan (5 — 2. These plots have been produced using the exact results given in the Appendix 
and not the MIA. 
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be obtained taking into account the bounds given by other processes. In Fig. 7 we plot 
distributions of those quantities for tan/3 = 2. 
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Figure 7: Distributions of X and Br{KL n^uP), Br{K'^ -k'^vD) for tan/3 = 2. 



As can be seen from Fig. 7, the allowed values for Br^Ki — > tt^i/u) and Br{K~^ — > 
TT'^i^p) can be significantly enhanced compared to the SM prediction, even by an order 
of magnitude. Such large values are higher than the 90%CL experimental bound [19], 
Br{K^ vr+z/z/) < 3.8 ■ 10"^'^, therefore this decay can already be used to constrain 
the MSSM parameter space. We will analyze this elsewhere. One should note that the 
obtained values for the branching ratios do not violate the Grossman-Nir (GN) [26] bound 
Br{KL — > 7r^iyiy)/Br{K~^ — > n'^uP) < 4.4, which can be regarded as a simple cross-check 
of the correctness of our numerical codes. 

In Fig. 8 we plot the allowed range of X after imposing the cut [19] Br{K~^ — > 
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Figure 8: Distributions of X for tau/S = 2 after imposing constraint Br{K~^ 7r~^i>i>) < 
3.8 ■ 10^^°. 



Ti^vv) < 3.8 • 10~^°. Even with this constraint, \X\ could be several times bigger than 
the SM value (2.4). Also, its phase can still vary almost freely, however is preferred to be 
in a broad range —160° <9x< 70°. Such a freedom leads to a possible enhancement of 
the ratio Br{KL — > 'k^vv) / Br{K^ n'^uP) even for large values of Br{K^ n^uP), as 
can be observed in Fig. 7. 

The plots of X in Figs. 7 and 8 display a conspicuous correlation between the phase 
and modulus. This can be understood from the allowed region in the complex plane for 
the chargino contribution to AjX shown in Fig. 4. The experimental constraints we apply 
constrain its imaginary part to be relatively small, especially when the real part becomes 
large. Taken together with the fact that Re At < (any 7) and that 7 is sufficiently 
constrained (cf. Fig. 3) to always imply ImAj > 0, this explains the shape of the allowed 
region in the complex X plane. 

Finally, in Fig. 9 we plot distributions of Br{KL n^uu) and Br{K^ ti^vv) 
for higher value of tan/3 = 20. Results are qualitatively similar to those obtained for 
tan/3 = 2, however for larger tan/3 (particularly tan/3 > 10) it is easier to generate 
parameters sets giving high branching ratios (or their ratio). 

To better illustrate the effects shown in the figures in this section, we give in Table 3 
specific examples of MSSM parameter choices leading to particularly interesting results. 
The quotation marks on some quantities in this table indicate that in evaluating these 
quantities we have set hadronic parameters to their central values as discussed in Sub- 
section 3.2. Varying these parameters, it is easy to reproduce the relevant experimental 
data within one standard deviation. 
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Parameter 


Example 1 


Example 2 


tan (3 


2 


20 


Ma 




2fi0 


1^ 


-375 


-344 




437 


928 


M2 


181 


750 


Msg 


308 


608 


Mi, 


138 


215 


Mt, 


279 


338 


Msi 


105 


884 


A 


-0.289 


-0.342 


7 


64° 


38° 



(2.18 - 5.02i) • 10-^ (7.57 - 0.871i) • 10"^ 
^ULR (-1.52 + 0.748i) • 10"^ 0.292 - 0.213i 



r23* 
'^ULR 


0.001 - 


- 0.604i 


0.239 - 


- 0.195i 


"1^ 1" 

\^k\ 


2.35 


•10-3 


2.10 


•10-3 


"AM/ 


3.15- 


10-13 


2.55- 


10-13 


"AM," 


1.03- 


10-11 


1.19 ■ 


10-11 


"5r(fi ^ X,7)" 


3.88 


■ 10-^ 


3.93 


• 10-^ 


Br{K+ TT+z/z/) 


1.78- 


10-10 


2.07- 


10-1° 




3.08- 


10-11 


4.34- 


10-1° 



Table 3: Examples of MSSM parameter points passing experimental constraints at the 
assumed accuracy (see discussion in section 3.2) and giving enhanced K'^ — > tt'^i/P and 
— > n^uP decay rates. Units are as in Table 1. See text for the meaning of quotation 
marks. 



23 



C; 

o; 

^/ 

O;' ... 



; o 

; -Q 

; Q- 
;-x 
; (D 
; _l ■ 
; O 

i O 



tan(S = 20! 



6 7 8 9 10 
Br(K"^TT"i/i/) -10-'° 



;^ 

*^ 4.5 

+ 4 

m 3.5 
\ 

;^ 3 

o 

^2.5 
is? 



m 



1.5 

1 

0.5 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


bound - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 
tan(S = 20 i 




. exp. 




'i'. ':*.•" 


; 1 

; 
; K 
; 
fen 
















»,?.*'.,;■- V. '. ■■ 
t>,ji ■ . 























6 7 8 9 10 



Figure 9: Distributions of Br{Ki^ — > t^^vv) and Br{K^ — > tt"*"!/!/) for tan/3 = 20. 



4.3 Comparison with earlier literature 

Colangelo and Isidori [30], extending and clarifying earlier work [38, 12], have used the fact 
that all contributions to the Z°-pcnguin must be at least of second order in the VEV by 
gauge invariance to trace the dominant SUSY contribution, at least in the case of large X, 
to the following term, which can be obtained as part of the chargino contribution (A. 14) 
to the Z'^-penguin: 

= -^^i^.^5.^^'*4''^'^'4'^^'^^*^r^rn"^*^^2(m^^,m^^,,m^J. (4.1) 

In the mass insertion approximation, this term contributes only starting from the second 

order, where it is proportional to the product Sjfu^Si^l^. 

Fig. 6 is consistent with this observation [30] , as is evidenced by the triangular excluded 
region in the bottom-right plot in that figure: |^| is only large when also \SijLjiSijlji\ is 
large. 

In fact, together with the large box contributions we reported in Subsect. 4.1 and Fig. 5, 
this plot suggests that even the box contributions, which were found to be subleading 
in [30] and earlier in [38], follow the same pattern, a point that has so far not been made 
in the literature. Indeed, there is no reason why the box diagram should not be competitive 
with the up-squark-chargino contribution to the Z°-penguin for light slepton mass. The 
chargino-box contribution (A. 7), ignoring slepton flavour violation and charged- lepton 
masses, contains a term 

W X T/ T/* ryRk* rySk rylm ryln* rylm ryln* 



Bl 



-dKLVRlVsj^U ^- ^- ^+ ^+ 



JIKL 2 

2 _2 2 2 



X mc^mc^DQ{mc^,mc^,mu^,mL^), (4.2) 

which is quite similar to (4.1). Consequently, it receives a similar second-order contri- 
bution in the mass insertion approximation, in contrast to the flnding in [30] that the 
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box diagrams contribute only at the level of three off-diagonal mass insertions. The dou- 
ble squark left-right flip is again the most important mixing parameter, even though no 
S'C/(2)-breaking powers of the VEV are required, unlike the case of the Z°-penguin. For 
the chargino contribution, that requirement therefore turns out to be a red herring, al- 
though it still explains why gluino and neutralino-penguin contributions are small. All 
large effects are proportional to the double up-squark left-right mixings simply because 
these are weakly constrained and at the same time avoid the double CKM suppression, 
both of which facts were correctly stated by the authors of [30] . 

More stringent bounds on the left-right mass insertions than those used in this work 
have been derived in [32] from a general MSSM RG analysis, however their results depend 
on the assumption that the soft-breaking terms are present already at the GUT scale. This 
is a model-dependent assumption, not necessarily true in the general MSSM. There are 
also potentially dangerous contributions from double mass insertions to AS — 2 processes 
that might lead to a violation of boTinds from — mixing. These are automatically 
taken into account in our computation by our exact diagonalization of the squark mass 
matrices. 

4.4 Extended scan over MSSM parameters 

The "adaptive scan" method [39] desribed in Section 3.3 allows for efficient exploration 
of really huge multi-dimensional parameter spaces, especially when the dependence of the 
analyzed results on most of those parameters is not very strong. Therefore, we tried to 
check how our results change if we get rid of virtually all assumptions usually used to 
relate MSSM parameters and treat them all as free independent quantities. We varied 
randomly the following quanties: 

• the angle 7 (real) 

• CP-odd Higgs mass Ma (real) 

• C/(l) gaugino mass Mi (complex) 

• SU (2) gaugino mass M2 (complex) 

• gluino mass rrig (real) 

• n parameter (complex) 

• diagonal left slepton mass m|,, common for all generations (real) 

• diagonal right slepton mass m|j, common for all generations (real) 

• 9 independent diagonal mass parameters in squark mass matrices, 3 parameters for 
each of left, up-right and down-right mass matrix (all real) 

• common sfermion LR mixing parameter A (real) 

• 3 independent LL mass insertions in squark mass matrices: S^^^, S^j^ (all com- 
plex) 
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• 6 independent RR mass insertions in squark mass matrices: Sj^^R, (^IxrR) ^^rri 
^URR^ (^URR, ^URR (all Complex) 

• 12 independent LR up- and down-squark mass insertions ^^^^^.i?) ^}dlrj ^^lri ^Idlri 

^DLRi ^IdLR-i ^IJLRi ^IJLR-i ^"uLRi ^IJLRi ^U^LRi ^u'lR (^U COmplex) 

Altogether this gives 63 real degrees of freedom, more than half of the free parameters 
of the completely unconstrained (but i?-parity-conserving) MSSM. For the SUSY mass 
parameters we use the same limits as given in Table 1, with the additional requirement 
I Ml I > 20 GeV; for new complex parameters we assume their phase to vary completely 
freely. Of course, for any parameter point we still apply all experimental and theoretical 
constraints listed in Section 3.2. Such an extensive scan was probably never reported 
before in the literature, so it is very interesting even to see if such a general MSSM 
version still retains any predictive power. The answer is positive - the results obtained 
for both analyzed decays, — > 7r°i/P and K'^ — > tt^uD, remain qualitatively similar to 
those already discussed for the "constrained" scan case! 

Plots equivalent to those shown in Fig. 3 confirm the assumption that the dependence 
on most new parameters is weak (i.e. distributions are fiat). Actually, in some cases like 
left stop and right sbottom mass parameters they become even flatter than in Fig. 3, as 
both discussed branching ratios do not depend directly on them (at least not strongly), 
and various experimental constraints can be satisfied varying additional parameters. 

The allowed ranges for X and both branching ratios are extended somewhat in the 
63-parameter scan, but not drastically, as illustrated in Fig. 10. In general, as could be 
expected, it is easier to generate a Ox in the full range, there are more points with large 
Br{K^ — > TT+i/P), Br{Ki^ n^uu) or their ratio, but in general all distributions are still 
well defined in shape, just a bit broader. Thus, the most important conclusion of this 
kind of analysis is the statement that even the almost fully general low-energy MSSM does 
not lose its predictive power and exploring it can still lead to reasonable and well defined 
results. Although it requires more effort in numerical computation, it also minimizes the 
possibility of overlooking some interesting scenarios which can be realized for particular 
SUSY parameter choices. 

The numerical analysis described in this section can also be used to put bounds on 
the allowed magnitude of the mass insertions in the squark sector, without resorting to 
the simplifying assumptions made in the traditional literature - absence of cancellations, 
neglect of interference of SM and SUSY contributions, the MIA. Results of such an analysis 
will be presented elsewhere. 

5 Conclusions 

In this paper we have analyzed the rare decays — > n'^uV and Ki, — > n^uV in a general 
MSSM with conserved R-parity. Our analysis goes beyond those present in the literature 
in that 

• we do not use the MIA but rather work in the mass eigenstate basis for all particles. 
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• we do not assume the dominance of certain contributions, as done frequently in 
the hterature, thus allowing for cancellations and interferences of different diagrams 
contributing to the decay rates in question, 

• we consider a scan of a very large space of parameters, 16 in the constrained scan 
and 63 in the extended scan, that to our knowledge has been presented here for the 
first time. 

We should emphasize that the physics scenario considered here differs from the simple 
new physics scenario discussed in [15, 16] in that the new physics contributions to X — > 
T^vv in the general MSSM do not only affect the function X as in [15, 16], but also modify 
the values of the CKM parameters. While the values of \Vus\. \Vub\ and \Vcb\^ measured 
in tree-level decays, are the same in both papers, the presence of relevant new phases in 
B^ g — B^ g mixings and Ek in the general MSSM can significantly alter the values of 7, (3 
and Rt that are used in the evaluation of the branching ratios for K — > ttuP. Consequently 
the comparison with the results presented in [16] has to be made with care. For instance 
the relation between f3x and 6x in (2.23) can be modified in the general MSSM through 
the change of the angle /?. While in [16] one has (3 = (23.5 ± 2.0)°, here one typically 
finds 12° < P < 27° but this modification can be compensated by a change of 9x so that 
Px remains unchanged. 

More important then is the difference between Rt used in [16] and Rt found here. As 
7 = (65 ± 7)° in [16] but 20° < 7 < 110° here (see Fig. 3), the corresponding values of 
Rt obtained from (2.19) with Rb = 0.37 are 0.86 < Rt < 0.95 and 0.66 < Rt < 1.18, 
respectively. 

While these differences in the CKM parameters play some role, they turn out to be 
sub dominant in comparsison with the differences between the values for X used here and 
in [16]: compare (1.9) and Fig. 10. 

Concerning individual contributions, in agreement with [12, 30, 32] we find that 
chargino diagrams are always strongly dominant, typically one order of magnitude larger 
than the neutralino contributions. Gluino exchanges are fully negligible. However, unlike 
these authors we find that in addition to chargino mediated Z^-penguins, chargino box 
diagrams can be important and even dominant for light charged slepton masses. Thus 
from the point of view of the decay rates Br{K^ ti^vv) and Br^Ki 7r°z/z/), the 
general MSSM is a new physics scenario with enhanced Z'^-penguins carrying a new com- 
plex phase, similarly to what has been considered in [16]. However, the presence of more 
free parameters than in the simple scenario in [16], allows for a broader range of values of 
|X| and 9x as we stressed above, and there are other manifestations of new physics such 
as in modified values of some CKM parameters and in box contributions to four-fermion 
operators. 

The answers to the questions posed at the beginning of our paper are as follows: 

• The phase Ox can be as large as found in [15, 16]. In fact as seen in Fig. 8 one finds 
typically 

-160° < < 50° (5.1) 

with a slightly increased range for the extended scan as seen in Fig. 10. Interestingly, 
there is a visible preference for negative values of Ox that seem to be required by 
the B — > t:K data within a simpler new physics scenario discussed in [15, 16]. 
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• However, among the allowed values of |X|, it is not easy to find simultaneously 
|X| ?a 2.2 and 9x ~ —85° as given in (1.9). Indeed for 9x ~ —85° the allowed value 
of \X\ in the constrained scan is typically lower than its SM value of 1.53 and only 
in the extended scan can it reach |X| ~ 2.0. Still, one cannot conclude from this 
finding that it is difficult to explain the B tiK data within the general MSSM, as 
the relation between the latter decays and K — > ttuu in the general MSSM is rather 
weak. A very recent study of S — > nK decays in the general MSSM can be found 
in [50]. 

• On the other hand for —50° < Ox < 50°, one can have \X\ as high as 7, implying 
that in this range of 9x very large departures from the SM expectations are possible. 
A similar situation is found for Ox — — (150± 10)° with \X\ reaching values as high 
as 7. This is illustrated in Figs. 8 and 10. 

• It should be emphasized that values of |X| significantly higher than 2, while excluded 
in the new physics scenario of [15, 16] through the data on S — > Xgl'^l" , are still 
allowed here because in the general MSSM there are no important correlations 
between B and K decays. 

• The pattern of Br{Ki^ ttVz/) ^ 3 • 10"^° and Br{K+ n+uu) ^ 8 ■ 10"^^ found 
in [15, 16], although not easy to achieve within the general MSSM in the constrained 
scan, as seen in Figs. 7 and 9, can be obtained in the extended scan as shown in 
Fig. 10. 

• More importantly, as seen in Figs. 7, 9 and 10, Br{K^ n^uu) in the ball park 
of the central experimental value in (1.4) can be naturally obtained and simulta- 
neously Br{Ki^ n^uu) can be even larger than Br{K'^ — > n'^uu). An example 
of supersymmetric parameters for which such pattern can be obtained is shown in 
table 3. 

• Comparing Figs. 7 and 8 we observe that the experimental bound Br{K^ — > n'^vD) 
in (1.4) has a significant impact on the maximal allowed values of |X| with essentially 
no impact on Ox- This is not surprising as n^uu is a CP-conserving decay- 
Clearly an improved upper bound on Br{Ki^ tt^vv) should have an important 
impact on the allowed values of Ox- 

• Finally, it is interesting to note that even in the general MSSM the angles /3 and 7 
of the unitarity triangle are rather constrained, 

12° < /3 < 27°, 20° < 7 < 110° . (5.2) 

The constraint from the measured asymmetry a^Xs (2.25) plays an important 
role here. 

In summary, within the general MSSM large departures from the SM expectations for 
K — > T^vv are still possible while satisfying all existing constraints. Br{K^ — > t^^vv) and 
Br{Ki, — > n^uP) can be both as large as few times 10~^^ with Br{Ki, — > n^uP) often larger 
than Br{K'^ — > tt'^uP) and close to its model independent upper bound. In particular 
the results of a phenomenological study of enhanced electroweak penguins in [16] can be 
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obtained. The supersymmetric effects thus turn out to be larger than found in [12, 32], 
where typically Br{K+ ti+uv) < 1.7 ■ 10"^° and Br{Ki^ vrVz/) < 1.2 ■ 10"^° have 
been found. With regard to [12], this is partly due to the fact that these authors use the 
single MIA, while we use exact mass matrix diagonalization. The second reason is the 
exploration of a much larger space of parameters that was possible in a reasonable time 
only by using specially designed Monte Carlo techniques. On the other hand, [32] took 
into account RG effects between the GUT and weak scales, as is suitable for a high SUSY 
breaking scale, which allowed them to constrain the double mass insertion responsible for 
the large effects reported in [30]. Working in the completely general MSSM, we did not 
impose the RG constraint, while our exact mass diagonalization incorporates the second 
(and higher) orders in the MIA, in effect reinstating these large contributions to the extent 
allowed by the constraints listed in Subsect. 3.2. 

In summary the large enhancements of Br{K^ ti^uu) and Br{Ki^ ti^uu) in the 
general MSSM found here are in the spirit of the findings of Colangelo and Isidori [30], 
except that we effectively incorporate not only the second but also higher order terms in 
the MIA, we explore much broader range of the space of supersymmetric parameters than 
done by these authors and we find the chargino box contributions more important than 
found in [30]. 
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Appendix Contributions to the function X 



In this appendix we collect the contributions to the functions and Xji within the 
general MSSM in the mass eigenstates basis for particles and sparticles. To compactify 
our notation, we abbreviate fermion- scalar vertices as iiV^Pi + V^Pr)-, with the actual 
values of and given by 
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where sw = sin 6*1^, cw = 0086*1^ and the definitions of the mixing matrices Z and the 
remaining vertices can be found in [49]. 

The contributions to Xl and Xji can be divided into two classes, generated by box-type 
and Z°-penguin diagrams: 

VsiV,*j[Xl,r]jikl= E {[BIr]jikl+[PIr]jiSkl}. (A.l) 

i=SM,h,c,n,g 

This formula corresponds to a generalization of (2.3) to 
3 

JIKL{djdi)v-A{'^Kl^L)v-A + [XR\jlKL{djdi)v+A{l^K'^L)v-A\ , 

K,L=1 

(A.2) 

where for K nuu decays one should use I = 1, J = 2 and K, L are neutrino fiavour 
indices (note that the subscript L on Xl refers to the chirality structure of the quark 
current and is no summation index). In general X^^ji quantities defined in (A.l) carry 
lepton fiavour indices. Only in the case when lepton fiavour number is conserved in 
the slepton sector and the left- and right-slepton fiavour diagonal mass parameters are 
identical for all three generations, the Xl^r are to a good approximation (neglecting small 
terms proportional to lepton Yukawa couplings) universal, as assumed in (2.3). Relation 
between X^^ji of (2.3) and of (A.2) is given by 

[Xl,r]21KL = Xlj^Skl- (A. 3) 

In the numerical analysis in our paper we assume mentioned above simple fiavour con- 
serving structure of the slepton sector and use the definition (A. 3). 
Defining further 

the non-vanishing box contributions read (summation over all indices other than J, /, K 
and L is understood) 

[Bl%iKL = -V,jV,^jAMxt) 5kl (A.5) 
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The non- vanishing penguin contributions read 
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The loop functions appearing in these Wilson coefficients are given by 
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(A. 17) 
(A.18) 
(A. 19) 

(A.20) 

(A.21) 
(A.22) 



(A.23) 
(A.24) 

(A.25) 



Infinite and //-dependent terms in C2 always cancels out in flavour ofi'-diagonal penguins 
after summation over squark, chargino and neutralino mixing matrices. We recall that 
fi{x) = Bq{x) and f2{x) = Co{x) in the notation of [7]. 
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